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Abstract. The one variable Bernstein-Szego theory for orthogonal polyno- 
mials on the real line is extended to a class of two variable measures. The 
polynomials orthonormal in the total degree ordering and the lexicographical 
ordering are constructed and their recurrence coefficients discussed. 



1. Introduction 
Let be a Borel measure supported on M^, such that 

f{x,y)dfi{x,y) < oo 

for every polynomial f{x,y). We would like to construct polynomials orthogonal 
with respect to this measure. Unlike one variable there are many natural orderings 
of the monomials a;"?/™ and each of these orderings gives a different set of orthogonal 
polynomials. Starting with [7], the preferred ordering is the total degree ordering 
and for polynomials with the same total degree the ordering is lexicographical, that 
is 

{k,l) <td (kiji) 

if 

k + I < ki + h OT {k + I = ki + h {k, I) <iex (kiji)), 

where lex means the lexicographical ordering (see below). In other words, we apply 
the Gram-Schmidt process to the polynomials ordered as follows {\,y,x,y'^,xy, 
x"^, . . .}. Let p!^ denote the polynomial of total degree n such that 

(1.1) K(^,y) = fc;'r'^^''y""''+ E ''^r^'y'^ 

{i,j)<td(r,n-r) 

with k^^~^ > satisfying 

p^-i^x^y-' dfi{x, y) — 0, 0<i + j<noTi + j = n and < z < r. 



(1.2) 
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Let Pn{x,y) be the {n + 1) dimensional vector with components the orthonormal 
polynomials of total degree n. The multiplications by x and y are given by three 
term recurrence relations 

(1.3) xPn = A^^r^-Pn+l + S^.„P« + A^,„_iP„-l 

(1.4) yPn = ^y.„P„+l + 5„,nP„ + 4,n-lP"-l- 

where A^^n, Ay^n are (n+1) x (n + 2) matrices such that rank(Aj;^„) — rank(Aj,^„) — 
n + 1, Bx^m By^n are symmetric (n + 1) x (n + 1) matrices and M* denotes the 
transpose of the matrix M. Notice that the block Jacobi matrices corresponding to 
multiplications by x and y will commute which amounts to certain commutativity 
relations between the matrices defined above, see [4] for more details. Since mono- 
mials with the same total degree are ordered lexicographically we see that Ay^n is 
"lower triangular" with positive entries in the i — 1 . . . ,n + 1 positions and 

A,j:_n is lower Hessenberg with positive entries in (i, i + 1), i — 1, . . . ,n + 1 entries. 

Recently an alternative way to approach two dimensional orthogonal polynomials 
was proposed in ^ by relating them to the theory of matrix valued orthogonal 
polynomials. This can be accomplished by using the lexicographical ordering, 

{k, I) <iex (^1, h) k < ki or (fc = ki and I < li), 

or the reverse lexicographical ordering 

{kj) <rcvlex (fcl,^l) {l,k) <iex (?l,fcl), 

to arrange the monomials. This naturally connects the theory of bivariate orthog- 
onal polynomials to doubly Hankel matrices, that is block Hankel matrices whose 
entries are Hankel matrices. 

For every nonnegative integer m we apply the Gram-Schmidt process to the basis 
of monomials ordered in lexicographical order i.e. 

• ■■,y"^,x,xy, . . . ,a;y™,a;^, . . .}, 

and define the orthonormal polynomials pl^ ^(a;, y) where < Z < to, by the equa- 
tions, 

(1.5) 

/ Pn m^*y"''^/^(^j y) = Oj < z < n and < j < m or i = ri and < j < I, 



and 

(1.6) Pi^{x,y)^k:il^,x-y'+ K^,r,.,i^y'- 

(ij)<lox(n,0 

With the convention fc^'^ ; > 0, the above equations uniquely specify Pn^m- Polyno- 
mials orthonormal with respect to d/i but using the reverse lexicographical ordering 
will be denoted by pj^ ^. They are uniquely determined by the above relations with 
the roles of n and to interchanged. Set 



(1.7) 



„o ■ 

Pn,m 
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The polynomials Vn,m may be obtained in an alternate manner as follows. We 
associate an (m + 1) x (to + 1) matrix valued measure dM™^^ (x) by taking 



(1.8) dM^+'ix) = JJ,i,y,...,yn'dK^,y)[i,y,---,yn, 

where the above integral is with respect to y. Let us denote by {P^}^=q the 
sequence of (to + 1) x (to + 1) matrix valued polynomials satisfying 

(1.9) P,™ {x) = ii:;j^„a:" + lower order terms, 

(1.10) / P,Tix)dM'^+'{PJP{x)Y = Sk,nlm+1 



with iC™„ a lower triangular matrix with strictly positive diagonal entries. The 
above conditions uniquely specify these left matrix valued orthogonal polynomials 
and it follows that 

(1.11) p„,™(x,y) = p,r(^)[i,y,---,yn'- 

From the relation between Pn,m and the matrix orthogonal polynomials we see 
that the following recurrence formulas hold 



where An^m and Bn^m are matrices of size (m + 1) x (to + 1) given by, 

(1.13) An^rn - {x¥n-l, 



(1.14) B„,™ - (xP„ 



m ; ^ n.m/ j 
m-j Pn , m ) . 



Here An^m is lower triangular with positive diagonal entries and Bn^m is symmetric. 
The analogous formula for P„^m(x, y) is 



(1.15) yP„,™ = A„,„+iP„.„+i + B 



n.ni^ n,m 



There has been much work done in studying bivariate orthogonal polynomials 
in the the total degree ordering and we refer to the books [B IH H] , and references 
therein. For Bernstein-Szego polynomials related to root systems see the newer 
article lOj. 

In this article we will consider a special class of "Bernstein-Szego" measures (see 
below for the definition) which are generalizations of the examples considered in [S] 
and construct orthogonal polynomials where the orderings used will be the total 
degree ordering or the lexicographical ordering. We will do this in order to inves- 
tigate and compare the properties of the polynomials constructed. The paper is 
organized as follows: In section 2 we briefly review the one variable Bernstein-Szego 
theory. We then consider two variable generalizations of this theory in section 3. 
Orthonormal polynomials in the total degree ordering as well as the lexicographical 
and reverse lexicographical orderings are considered. In section 4 we investigate 
the properties of the coefficients in the recurrence formulas satisfied by these poly- 
nomials and in section 5 we consider several examples including the ones discussed 
in[S]. 
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2. One variable Bernstein-Szego measures 

Consider a measure supported on [—1, 1] absolutely continuous with respect to 
Lebesgue measure of the form d^{x) — w{x)dx with w{x) = (1 — x'^)^/'^ / p{x) where 
p{x) is a polynomial of degree I positive on [—1, 1]. Such measures are in the class of 
Bernstein-Szego measures which have many nice properties. We call a polynomial 
h{z) stable if h{z) ^ 0, \z\ < 1. We have 

Theorem 2.1. Let w{x) = (1 — x^Y^^ / p{x) where p(x) is a polynomial of exact 
degree N and is positive on [—1,1]. Let p{x) = |/i(e*^)p where h{z) is a stable 
polynomial of exact degree N in z such that h(0) > 0. Then 



2X1/2 



Pn{x)^i^-j {sin9)-^'^{e"^''+^>'h{e'')), x ^ cos9, N < 2{n + 1), 

is an orthonormal polynomial of degree n associated with w. 

The above Theorem shows that the orthogonal polynomials associated with 
Bernstein-Szego weights can for large enough n be written as a linear combina- 
tions of Chebyshev polynomials (for instance [5]) . A proof of the above Theorem 
in the case where the coefficients of h are real is given in Lemma [3T] For real coef- 
ficients the above formula holds for = 2(n+ 1) if the right hand side is multiplied 
by (1 — hN/ho)~^^'^ , where /lo = h{0) and /ijv is the coefficient of z^ in h{z) We 
now show how to construct orthogonal polynomials of degree < [-^^y^] — 1. 

Proposition 2.2. Assume that h[z) = 1-f hiz + ■ ■ ■ + hjsjz^ is a stable polynomial 
of degree N with real coefficients. Let qu be defined by 

N 

(2.1) gfe(a;) =^/i,C/fc_,(x), /c > 0, 

4=0 

where Un{x) = —U-n-2{x) for n < 0. For < k < [ "'^~^ ] — 1, there are constants 
h'l. such that 

qk{x) qk{x) + h'f^^^qk+i{x) H h ^^_fe+29w-fc-2(a:) 

is a polynomial of degree k and orthogonal to every polynomial of degree less than 
k with respect to 

\f~\~^'x^ dx 

= \U( \\2 ' x = (z + l/z)/2, 

where ^^_j,^2 ^ ^N-k+i ~ ^-/v-i — hjqhi, and the other h'j can be deduced 
inductively from hj . 



Proof. Lemma [3TT] shows that qk is orthogonal to all polynomials of degree at most 
k — 1. Thus, it follows readily that for < k < [-^^y^] — 1, qk is orthogonal to all 
polynomials of degree at most k — 1. We now prove that we can choose constants 
h'j such that qk is of degree k. Using C/„(x) = —U-n-2{x) for n < 0, we can write 
(2.2) 

Af-fc-2 k \N-2' 

qk{x) = - ^ hk+i+2Ui{x) +^{hk-i - hk+i+2)Ui{x), k< — — -1, 

i=k+l i=0 
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in which the first sum contains terms that have degree > k. In particular, degqk < 
N — k — 2 for k < f'^^j^] — 1- Furthermore, for fc + 1 < < — 1 we can write, 
by (|2.ip and the fact that Un{x) ~ ~U-n-2{x) for n < 0, 

3 k N-j-2 

Iji^) ^ X! hj^iUi{x) +^hj_iUi{x) - ^ hj+i+2Ui{x), 

i=k+l 1=0 i=0 

where again the first sum contains terms that have degree > A: in x. Since ho{y) — 1, 
it follows readily that 

N-k-3 

qk{x) + hNqN-k-2{x) = - ^ {hk+z+2 - hNhN-k-2-i) Ui{x) H , 

i=k+l 

where only terms whose degree > A; + 1 are given explicitly in the right hand side. 
Thus the right hand side of the above expression has degree N — k — 3. Continuing, 
we add (/ijv-i — ^Jv/ii)9iV-fc-3 to eliminate UN-k-3- Proceeding in this way, we 
keep adding terms until the right hand side contains only terms of degree < k. This 
proves that qk is indeed a polynomial of degree k. □ 

Example 2.3. The first non-trivial case of Theorem 12.11 and Proposition 12.21 ap- 
pears when N — 5 for which f-^^^y^] — 2 and we are missing the orthogonal poly- 
nomial of first order. In this case, the formula (|2.ip shows that 

qi{x) = ~h5U2{x) + {ho - hi)Ui{x) -f {hi - hj,)Uo{x) 
q2{x) = hoU2{x) + {hi - h5)Ui{x) + {h2 - hi)Uo{x). 

Our qi is given by 

qi{x) = qi{x)+hzq2{x) = [(/ii-/i3)+/i5(/i2-/i4)]?^i(a;) + [(/io-'j4)+/i5(/ii-^5)]C^o(a;), 
indeed a polynomial of degree 1. 

Remark 2.4. The above proposition shows how a complete orthogonal basis can 
be derived for the Bernstein-Szego weight function. A similar approach will be 
used in the construction of bivariate polynomials orthogonal in the lexicographical 
ordering. However, this procedure will not help us with the construction of bivariate 
polynomials orthogonal in total degree ordering. Indeed, when hk are polynomials 
of y, the degree of qk as a polynomial of two variables could have total degree 
greater than k. 

We will also investigate the recurrence coefficients associated with Bernstein- 
Szego weights. Let ^.{x) be a positive measure supported on the real line with an 
infinite number of points of increase and {pm\ be the set of polynomials of degree 
m orthonormal with respect to each with positive leading coefficient. Then it is 
well known that they satisfy the recurrence relation 

XPm{x) = ara+lPm+l{x) + bmPni{x) + 

If the measure /i is of the form considered in Theorem 12.11 then the formula for 
Pm above shows that for I < 2m they satisfy the same recurrence formula as the 
Chebyshev polynomials i.e a„i = 1/2 and bm = 0. For extensions of this result see 
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3. A CLASS OF Two Variable Bernstein-Szego polynomials. 

Here we study two variable polynomials orthogonal with respect to a Bernstein- 
Szego weight of the form 



pix, y) 

where p is a polynomial in x and t/, positive for — 1 < < 1 and — 1 < y < 1. 
Lemma 3.1. Let N be fixed positive integer. For any i = 0, . . . , N let hi{y) be 

2 I 2 



polynomials in y with real coefficients of degree at most ^^'^ f^o{y) = 1; 



such that 

N 

(3.1) h{z,y) = Y,hi{y)z\ 

i=0 

is a stable polynomial in z for all — 1 < y < 1, i.e. h{z,y) ^ for any \z\ < 1. 
Define 

N 

(3.2) qk{x, y) = ^ hi{y)Uk-i{x), 

i=0 

where Un{x) is the n-th Chebyshev polynomial of the second kind. Here, if n < 
the Chebyshev polynomial is understood as Un{x) = —U-n-2{x). 

Then qk{x,y) is a polynomial in two variables, which is orthogonal to every 
polynomial in x of degree less than k with respect to 



dfJ-yix) = ^^^^ y)|2 ^^' x=iz + l/z)/2. 
Moreover for k > <lk{x,y) is a polynomial of total degree k and 



L 



\l{x,y)d,iy{x) = {^ ^/fc^rr'^a'l 



Proof. We begin by showing that qk is of total degree k when k > \^^-^^ • Note 
that |/iAr(y)| < 1 follows from the stability of h. If fc > A^, then qk is of degree k 
follows obviously from — |^ — «| < i. Let now k < N — 1. Then using the fact 
that U-j-2{x) = —Uj{x) for j < 0, we have 

k N 
(3.3) qk{x,y) = ^hi{y)Uk-i{x) - ^ hi{y)Ui-k-2{x). 

i=0 i=k+2 

It is easy to see that §-\^-i\+i-k-2<N-k-2<kask> \^]. Hence, 
Qk is of degree k. 

To establish the claimed orthogonality properties of qk we follow Szego. Since 
Un{x) = - — jEfji — ) the polynomials qk{x, y) can be written in the following form 

1 ^ 

= [z''+'h{l/z,y)-z-''-'h{z,y)] , 

and this identity holds for every k. 
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Take j < k and integrate qk{x, y)Uj{x). Substituting x = cos 6 we get 

io {z^l/z)^h{z,y)hil/z,y) ^ ' 

^ - '-"'^ odo - r -:^'^!-^:;-r';^ sin^ odo. 



iz-l/zmz,y) ''''''' y„ {z-l/znil/z,y) 

If ^? is replaced by —9 in the second integral, the two integrals can be combined into 
one integral over [— 7r,7r]. This can be rewritten as a contour integral over \z\ = 1, 
with z = e'^ as follows 

Since h{z, y) is a stable polynomial the last integral vanishes for fc > j by the 
residue theorem. 

Finally, let us see that for k > the norm of Qk does not depend upon y. 

First let TV < 2fc + 2. Using ((H^ we obtain 

qk{x,y)qk{x,y)— rr^da; = / qk{x,y)hQ{y)Uk{x)--- -j^dx 



1 



Hz^yW 7-1 '"^ '\Hz,y)\ 
thoiy) f zJ>'+^ - z-^ -ihojy) f z'^ 



4 J\z\=i h{z,y) 4 Ji^i^i 

" 2/i(0,z/) ^ 2- 

For the case = 2fc + 2, p.3p shows that qk{x,y) = {ho — hN)Uk{x) + so 
that a simple modification of the previous computation gives that the norm of qk 
is ^il-hNiy)/hoiO,y)) = ^{l-hN). □ 

Remark 3.2. The condition deg hi{y) < y ~\y ~ *l Lemma [Oi ls sufficient for 
the class of Bernstein-Szego measures considered in the present paper. However, 
the statement can be easily extended to more general weights. For instance, if 
deg hi < i for every i = 0, 1, . . . , A^, then qfc(x, y) is a polynomial of total degree k 
in X and y for fc > iV — 1. 

Let n„ denote the space of polynomials of total degree at most n and n„ „ = 
span{x'y^,0 < i < 7t.,0 < j < m}. Using the above Lemma, we can construct 
polynomials of total degree n orthogonal to all polynomials in n„„i. 

Theorem 3.3. Consider the two variable measure 



3.5) dfi(x,y)^— — — dxdy, 

TT^ \h{z,y)\^ 

where h{z,y) is as in ()3.1|1 . Then, for ^ k < n, the polynomial 

Pn{x,y) = qk{x,y)Un-k{y) 

is orthogonal to all polynomials of degree less than (fc, n ~ k). 
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Proof. From Lemma l3.ll it follows that qk{x,y)Un^kiy), ["^^2^1 < ^ < is of 
degree (fc,n — k) in n„. The result follows directly from previous lemma and 
the orthogonality of Chebyshev polynomials. To see this consider the product 
Ui{x)Uj{y) \ii < k,i+j < n and qi{x,y)Uj{y) if i > k,i + j < n. Then, integrating 
first with respect to x and using the previous lemma yields, 

1 

Pn{x,y)Ui{x)Uj{y)dn{x,y) = 0, 

-1 J-i 
a i < k and also 

fi pi 



Pnix,y)qiix,y)Uj{y)dfi{x,y) =0, 

-1 J-i 

for i > k. Finally, for i — k we perform the integration with respect to a;, and 
since the norm of qk{x,y) does not depend on y, the orthogonality of Chebyshev 
polynomials shows that for k + j < n, 

2 



Pni^^,y)qk{x,y)Uj{y)dfi{x,y) ^ - / Un-k(y)Uj(y)\/l - y'^ dy ^ S^-k,: 



The result follows since if x^y^ £ n„ is a monomial of degree less than (fc,n — k) 
then 

x^y^ e spa.n{Ui{x)Uj{y),{iJ) <td (fc,n-fc), i < k} 

e span{(j,;(x,?/);7j(2/), (i, j) <td (fc,7i - k), i > k}. 

□ 

In a similar way, if 

M 

h{x, w) ~ hj {x)w-' , 

3=0 



is a stable polynomial in w for any — 1 < a; < 1 where hj(x) are polynomials with 

'jV-^; — "2 I "2 



real coefficients such that the deg hj{x) < ■¥■ — — j| and ho{x) = 1, we have: 



Lemma 3.4. For any I such that \ ^'^^ ^ ] < I the polynomial 

M 

qi{x,y) = ^hj{x)Uk^j{y), 
is orthogonal to any polynomial in y of degree less than I with respect to 



duxiy) = ^ ^ -dy. 
^ ' \h{x,w)\^ 

Moreover the norm of qi is independent of x. 

Analogous results hold if we order the monomials in the lexicographical or reverse 
lexicographical orderings. Let T\n,m be given as above and order the monomials 
according to the lexicographical ordering. Let k = maxdegree(/ii(y)) and qr be 
given by equation (|3.2p . Then we find 
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Lemma 3.5. For ["^^j^l < r <n and < k < m — k, 

Pr.m = qr{x,y)Uk{y) 

is of norm one, of lexicographical degree {r,k), and is orthogonal with respect to 
the measure (j3.5p to all the monomials in Tln,m of lexicographical degree less than 
(r, k). When N is even and r = N/2 — 1 is given by the above formula up to 

a multiple. 

Proof. From the definition of qr{x,y) and the constraints on the degree of hi{y) 
we see that qr{x,y)Uk{y) G for k < m — k and is of lexicographical degree 

j-Tyk^ The orthogonality of qr{x,y)Uk{y) now follows from Lemma [3.11 and the 
orthogonality of Uk{y)\/T—y^ to y^ , < j < k. □ 

The above Lemma shows the remarkable fact that for the weights we are con- 
sidering if we increase m, ^ does not need to be recomputed. 
Likewise with 7 = maxdegree(/ij(a;)) we have 

Lemma 3.6. Let 

pLt = qt{x,y)Ui{x). 

Then for \ '^'^^ ^ ] < t < m and 0<l<n — ^,pl^^ is of norm one, of reverse 
lexicographical degree (l,t) and is orthogonal with respect to the measure p.Sp to 
all the monomials in Tln.m of reverse lexicographical degree less than {l,t). When 
M is even and t — M/2 — 1 pj^ j is given by the above formula up to a multiple. 

At this point without any other assumptions on a in both orderings there are an 
infinite number of missing polynomials. More precisely in the total degree ordering 
we are unable to compute for k < while in the lexicographical ordering 

Prm is not given for r < f-^^j^l or k > m — k. We now consider a class of weights 
where we can compute „ for k > m — k . 

Let 

3.6) (T{x,y}^— — = — — , 

\Hz,y)\ \h{x,w)\^ 

where y) — uj{z, u))uj{z, 1/w), with 

N 

(3.7) a'(z, t(j) = JJ(1 + a^zw), 

1=1 

x = i(z + 1/z) and y — ^{w + 1/w). Here we assume < \ai\ < 1, a.i <E M. 
Thus Lu{z,w) is a stable polynomial, i.e. uj{z,w) 7^ 0, \z\ < 1, \w\ < 1. Since for 



z 



1 = |w| we have \h{z, y)\ — Y[i=i{l+o,iZw){l+aiz/w){l+aiw/ z){l+ai/{zw)) 
we see that 

(3.8) h{x,w) = uj{z,w)uj{l/ z,w). 

We begin with 

Lemma 3.7. The polynomial w^uj{z, 1/w) is homogeneous in z and w of degree 
N . With h{z, y) = X)i=o ^i{y)z^ then deg hi < N — \N — i\ and h{z, y) is stable for 
— l<2/^l)|-z|<l- The same properties hold for h with the roles of z and w and 
X and y interchanged respectively. 
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Proof. The homogeneity of w^a>(z,l/w) follows from the definition of uj{z,w). 
Since (1 + aizw) is stable gi{z) = (1 + aizw){\ + aiz/w) = (1 + 2aiyz + afz"^) 
is stable for —1 < 2/ < 1 and \z\ < 1. Furthermore with gi{z) ~ X]j=o 3«J (j')'^"' 
we see that deg gij — 1 — |1 — j|. The degree properties of hi{y) now follow since 
h{z, y) — ni=i 9i{^-i y)- analogous argument holds for h. □ 

Lemmas l3.5l and l3.6l allow us to compute pjf „j and ^ for r and t sufficiently large 
and k < m — N and I < n ~ N respectively. We now embark on the construction 
of „j and J for m — N < k < m and n — N < I < n respectively. In order 
to accomplish this, let K = span{Un{x)Um{y) , n > 0, m > 0} i.e. the space of 
all polynomials in x and y, and K = span{z~"ii;~'", n > 0,to > 0}. We define a 
linear map T : K ^ K T{Ui{x)Uj{y)) — z~^w^^ where y = l/2{w + 1/w) and 
X ~ l/2{z + 1/ z). It is not difficult to check that T is a linear one to one mapping 
from K onto K. We begin with the simple 

Lemma 3.8. For fixed n and m let p{x) be a polynomial of degree at most n in x 
and q{y) be a polynomial of degree at most m in y. Then T (jp{x)q{y)Un{x)Um{y)) = 

Proof. Suppose deg p{x) — p and deg q{y) = q. Then from the recurrence for- 
mula for Chebyshev polynomials i.e., Uk+i{t) + Uk-i{t) = 2tUk{t) we see that 
p{x)Un{x) = Y^^j=-p'^j^n-j- From the restrictions on the degree ol p it follows 
that the indexes of the Chebyshev polynomials in the above sum are positive so T 

can be applied to obtain T{p{x)Un{x)) = X)j=-p ^i^"""'^"' = z~'^p ^^-^^^^ Since 
the same is true for q(y)Um{y) the result follows from the linearity of T. □ 

The purpose of the map T is that it conveniently keeps track of the highest 
powers of x and y in terms of z and w. 
We now prove. 

Theorem 3.9. Let p'^^{x,y), m — N < k < m, with m,r > 2N be of norm 
one, of lexicographical degree (r, k) and orthogonal to all monomials in Hn.m of 
lexicographical degree less than {r,k). Then there exist constants am-k,j, bm-k,j, 
j — 0, . . . , k ~ (m — N) ~ 1 such that 

k-{m-N)-l 

{x,y)Uk-^j{y) + j{x,y)Ur+ 

with arn-k,Q ^ 0. 

Proof. From Lemma 13.11 we see that q,„+i+j(a;, is orthogonal to all monomials 
y*, < i < m and qr+j{x, y), j = 1, . . . , k — {m — N) — 1 is orthogonal to all x^ , I — 
0, . . . , r. Furthermore qr{x, y)Uk{y) is orthogonal to x^y^ for Q<i<r,Q<j<m 
and for i = r, < j < k. Thus the right hand side of the above equation satisfies 
the required orthogonality conditions. We now show that it is a linear combination 
of the monomials of degree less than or equal to (r, k) in lin.m ordered lexico- 
graphically. Set Sr,ra-N+1 = qr{x,y)Um-N+liy), Sr^N,m+l = qm+l{x , y)Ur-N (x) , 

Sr,m-N+i = T{qr{x,y)Ura-N+i{y)) a.ndSr-N,m+i = T {q,yi+i{x , y)Ur- n{x)) . Then 



ON A TWO VARIABLE CLASS OF BERNSTEIN-SZEGO MEASURES 



11 



from Lemmas 13.71 and 13.81 we find that 



■2N / , 1 / 

C -r -(»ri-JV+l) . W-tl/W 



2 

i=0 ^ 

which contains i/;^'-™^"'^' . Tlierefore we must eliminate the w*^ term in ll!(z,w)w^uj{z, l/w). 
From Lemma 15771 we have that this is associated with z^ . Now 

2^ 



Sr~N.m+l = W ™ ^ ^ hi 



z + 1/z 



w'z-'^+^ 



, 2 

W-"'-^Z-'-Lu{z,w)z'^Uj{l/z,w), 

where equation (|3.8p has been used to obtain the last equality. From the definition of 
uj{z, w) we see that the coefRcient of z^w^ in z^ lo{1/ z, w) is 1 and that z^lu{1/z, w) 
and w^uj{z, l/w) are homogeneous polynomials of degree in z and w. Write 



and 



Z^L0il/z,w)^Z^ + ---+jlfW^. 

Then Sr,m-N+i — kj^Sr-N.m+i docs not contain if kj^ = jj^. Thus the 

polynomial = Sr^m-N+i^ kj^^Sr-N.m+i is of lexicographical degree {r,m — 

N+1) and is orthogonal to all polynomials of degree less than {r,m~ N +1). Note 
that 

Sr.m-N+l - kljSr-N,rn+l = Z^"^ W^"^ Ujiz , w)uj^ {z , w) 

where (^^{z, w) is a polynomial homogeneous of degree iV — 1 in z and w. 
To continue, set 

^r,m ~ Sr,m — N+1 ^ k^ Sr-N,m+1 

and 

From the above discussion we see that T{S^„J = z~^w'"^ll!{z,w)ll!^{z,w) where 
uj^{z,w) is a polynomial homogeneous in z and w of degree N — 1, i.e., 

2 = 

with 7g 7^ since is of degree {r,m - N + 1). Therefore T{Sl„,+^) = 

z^^w^"^^^uj{z,w)uj^{z,w). Likewise because of the relationship between Sr,m and 
Sr,m we find T^CiS^ m+i) = z~'^w~™~^uj{z, w)oj^(z, w) where uj^{z,w) = z^~^w'^~^uj^{l/ z, 1/ 

If we choose k%_^ = thenT{S^„^_^^-k%_^S^„^^^) does not contain . 

Hence 

~m-JV+2 _ cl _ , 2 cl 

is of lexicographical degree {r,m — N + 2) and is orthogonal to all monomials in 
n„,m of lexicographical degree less than {r,m — N + 2). 
Set 

c2 _ cl _ j,2 cl 

and 

c2 _ cl _„ 1,2 cl 

'-'r,m ~ ^r+l,m '^N-l'-'r+l,m- 
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The above discussion shows that 

with uj'^{z,w) a polynomial homogeneous in z and w of degree N — 2 with the 
coefRcient of being nonzero. Likewise T{S^j-,J — z'^^w~"^uj{z,'w)uj'^{z,w) 

with w) = z^~'^w^~'^uj'^{l/z, 1/w). The construction of an orthogonal set of 
polynomials of the correct degree now follows by induction, that is, 



k-{m-N)-l 

E 



To complete the proof we show that the norm of p^m depends only on m — fc. 
From the definition of q and q and Lemma |3. II we see that 

Pr.mi^i v) = Crn~kUr{x)Uk{y) + tcrms of lowcr lex degree. 
From the orthogonality properties of q and q we find, 

llPr.mlP = C,n-k{Pr.m,Ur{x)Uk{y)) = C,n-ka,n-k,Q {qr{x , y)Uk{y) , Ur{x)Uk{y)) ■ 

By orthogonality Ur{x) may be replaced by qr{x, y) showing that mll^ — Cm-kam-kfi 
which proves the result. □ 

A similar discussion shows 

Theorem 3.10. Letplf{x,y), n-N <t < n, n,t > 2N be of reverse lexicographi- 
cal degree orthogonal to all monomials inlln,m of reverse lexicographical degree 
less than {I, t). Then there exists constants at-ij, bt-i,j , J = 0, . . . , Z — (n — N) — 1 
such that 

i-(n-Af)-l 

Pn,t{^^v) E {at-Ljqt+A^^v)Ui~i{^) + ~^t-Ujqn+i+j{x,y)Ut+i-n~i-j{x)^ . 

3=0 

The above discussion gives, 

Theorem 3.11. Let a{x,y) be given as in equations p.6p and p.7p . Let Tln,m — 
span{a;*y-', <i<n;0<j<m} with n > 2N and m > 2N. Then for 2N < r < 
n, < k < m 

qr{x,y)Uk{y) if k < m - N 

I]j=o'"~^''~"^ {0'm-k,jqr+j{x,y)Uk^j{y) + bra-k,jqm+l+j{x,y)Ur+k-'m-l-j{x)} 

if k > m — N. 
Likewise for 2N < t < m 

PnA^^y) = 
' qt{x,y)Ui{x) if l<n~N 

Ej=o"~^''^^ [at-Ljqt+jix,y)Ui-j{x) + bt-Ljqn+i+j{x,y)Ut+i-n-i-jiy)^ 
if l>n~ N. 
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4. Recurrence Coefficients 



We now examine the consequences for the recurrence coefficients if we assume 
that the weights are of the form discussed above. From the one variable theory one 
might expect that the recurrence coefficients have some simple structure. We begin 
with the total degree ordering. 

Theorem 4.1. Consider the weight given by equation ()3.5p . For \ ^2^ '\ < n, 



(4.1) 



-^y,n ■ — 



1/2 







1/2 0_ 



where Cy^„ is a \ ^2^ ^ ^ [ "^^2^ 1 ^ower triangular matrix with positive diagonal 
elements 



B,, 



Dy.n 









(4.2y -^y,n 



where Dy n is a symmetric f ^^*^ ] x [ ^^'^ ~\ matrix, 

'C 



(4.3) 







l/2/„_piv_^^_„_^iv+i| 



where Cx,n is an \ ^2^ ~\ x ['^^j-'-] lower Hessenberg matrix with positive entries in 
the upper diagonal and 



(4.4) 



Dx,n 





where D^.n is a symmetric [^] x [-y] matrix . Here ^ [ "+'- ] '^'^ 

{n — \ ^2^ ^) X (n — f-^^^l) matrix with ones on the {i,i + 1) entries and zeros 
everywhere else. 

Proof. If we examine the matrix elements in Ay,n we see that for [ ''^^"'" 1 ^ h 



[^i/,n]j,j 



yqi{x, y)Un-i{y)pi+i{x, y)d^{x, y). 



The recurrence formula for the Chebyshev polynomials shows that yqi{x, y)Un-i{y) - 
i(Ki-i(^'y) +Pn+i(^'y)) which implies equation (|4.ip . Equation (|4.2p follows in 
a similar manner. To show (|4.3p we write for f-^^^-'-] < i 



[A, 



xqt{x, y)Un-t{y)pi+i{x, y)dii{x, y). 



From the definition of we see that a;gi — ^(gi-i+gi+i). SotYidX xqi{x^y)Un-i{y) = 
^(p^l\(a;, y) + P^n+ii^^y)) which gives (|4.3p . Equation (|4.4p follows in a similar 



manner. 

For the polynomials ordered lexicographically we have, 



□ 
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Theorem 4.2. Consider the weight given by equation ()3.5p . In the lexicographical 
ordering we have for \ ^2^ ~\ < n 



(4.5) 



A — 



1/2/m-re+l 
Cn.T 



where Cn.m is a k x k lower triangular matrix with positive diagonal entries and 
(4.6) B„,™ 






with Dn.rn « Symmetric k x k matrix. For weights of the form p.6p - p.7p and 
2N < n,m we have 



(4.7) 

and 

(4.8) 



An,m = 1/2/m+i and B„ 



Ai,7n = l/2/„+i and B„ 







Proof. The matrix elements in An^m are given by 



[An, 



Thus equation (14. 5p follows from Lemma [3751 If the weight is of the form (|3.6p - (|3.7p 
and 2N < n,m then the first part of equation (14. 7p follows from Theorem l3.9l The 
remaining claims follow in a similar manner. 

□ 

5. Examples 

In this section we consider several examples including the ones discussed in ^ . 

5.1. Example 1. Here we consider the polynomials orthogonal with respect to the 
probability measure on [—1, 1]^ 

4 



(5.1) 
where 
(5.2) 



dfi 



w{x,y) 



tw{x, y)V^ - x'^\/i - y'^dx dy, 



(l-a2) 



4a2(a;2 + y^) - 4a(l + a^)xy + (1 - 0^)2 ' 

and a is real and \a\ < 1. This measure is of the form equations (|3.6p and (j3.7p 
with iV = 1 and h{z, y) = 1 ~ 2ayz + a^z^. 

From Lemma l3.1l we know that the polynomials qo{x,y) — 1 and 

qk{x,y) = . ^ „ {Uk{x) ~ 2ayUk-i{x) + a^Uk-2{x)) 
V 1 — a 

for fc > 1 are orthonormal with respect to the measure 



dny{x) = -- 



{1-0^)^/1^ 



TT 4a2(a;2 + y2) _ 4^(1 _^ ^2)3.^ ^. „ ^2)2 

Applying Theorem 13.31 shows that for fc > 
(5.3) 



dx. 



Pn{x,y) = qk{x,y)Un-k{y)- 
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This gives the recurrence coefficients [5] 



A — - 



a \/l — 



1 



and i?a;,n, By^n are (n + 1) x (n + 1) zero matrices. 

For the polynomials obtained via the lexicographical ordering Theorem 13.111 
shows that for n, m > 2, 



(5.4) 



Pn 



qk{x,y)Uk{y) for /c < TO 

aiqn+i{x,y)Um{y) + biq,n+iix,y)Un^i{x) for k ^ m. 



5.2. Example 2. Consider the weight w{x,y) associated with 

h{z,y) = {l-2bz)(l-2ayz + a'^z'^), |6| < 1/2, |a| < 1. 
Lemma |3 . 1 1 shows that for fc > 1 the polynomials 

(5.5) qk{x,y) = Uk{x) - 2{ay + b)Uk-i{x) + a{Aby + a)Uk-2{x) - 2a%Uk^^{x) 
are orthogonal with respect to the probability measure 



dfiy{x) 



2 w(x, y)y/l ~ x-^ 



TT (1 - 46x + 462) 
on [—1, 1], where w{x, y) is defined in (|5.2p . It is easy to see that 



(5.6) 



dlly{x) = — 

1 



1 



\b\ < 1/2, 



iaby + Aa^h^ ' 

which is another Bernstein-Szego weight whose orthogonal polynomials are easily 
seen to be 

'Un{y)-2abUn-i{y), if |2a6| > 1 
2abUn{y)-Un-i{y), if |2afe| < 1. 
Hence, applying Theorem 13.31 and defining qo{x,y) = 1 shows that 



Vn{y) 



(5.7) 



Pnix^y) 



qk{x,y)Un~k{y), i <k <n 



give a complete set of orthogonal polynomials in the total degree ordering. It follows 
that the recurrence coefficients in the total degree ordering are as suggested in j5], 
i.e. Ax^n and Ay^n are the same as in the previous example. For i?x,n we have 



Bx,o — [b], Bx^i — b 
and for n > 2, Bx,n is the block matrix 

Bx'n — 



— aVl — a? 



Bx.i 




-a\/l — 

^2 
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The matrices By,n are identically equal to zero for n > 1 and Byfi = [6a]. 

5.3. Example 3. For completeness we include also the formulas for the second 
example considered in [5] . The measure in this case is given by 

2zo 



(5.8) " , 

^ ' 2(1 — z ) 

H —w{x,y)6{x — xo)\/l — x'^\J\ — y^dxdy, 

TT 

where again w{x,y) is defined in (j5.2p . Here zq = is real with magnitude less 
than one, xq — j{zq + Zq^)^ the first part of the measure is on [—1, 1]^, the second 
is on M X [—1,1] and w{x, y) is given by (|5.2p . 

If |5| < 1/2 then w(cc, y)\/l — x"^ is exactly the Bernstein-Szego weight 

discussed in Example 2. When |6| > 1/2, /i is no longer stable and a modification 
of the calculation in Lemma 13.11 shows that the polynomials given in (j5.5p are 
orthonormal with respect to the measure 

d^y{x) — — — rw{x, y)y/T~x^dx + (1 — Zq)w{x, y)S{x — xo)y/l — x'^dx. 

7r(a;o - x) 

A residue calculation gives, 

(5.9) I dtiy{x) ^ 



J y^^' l + 4a262_4a5y' 

which is the same as (j5.6p . Hence, it follows that the polynomials given in (|5.7p are 
also orthogonal with respect to d^{x, y) in (|5.8p . 

The polynomials and the recurrence coefficients in the lexicographical ordering 
are given in [5] using a Darboux transformation connecting the two examples. 

5.4. Example 4. In this example we consider 

h{z,y) = (1 - 2aiyz + alz'^){l - 2a2yz + a^z"^), |ai|, |a2| < 1. 
with corresponding measure, 
4 



dfi{x,y) = —w{ai;x,y)w{a2;x,y)\/l - x'^^/l - y'^dxdy, 

in which 'w{a;x,y) is defined as in (15. 2p . Lemma [3.11 shows that for fc > 1 the 
polynomials 

qk(x,y) =Uk{x) - 2aia2(ai + a2)Uk-i{x) + {af + aj + 4ai02y^)J7fe_2(a;) 

- 2aia2(ai + a2)Uk-3ix) + a^b^Uk-iix), 

where U^m~2(x) ~ —Um{x), are orthogonal with respect to the measure diJLy{x). 
On the other hand, a simple computation via residues shows that 

da (X) = ^ + "^^2 

-1^"^ ^ (l-aia2)((l + aia2)2-4aia22/2)- 

This is again a Bernstein-Szego weight, whose orthogonal polynomials are given by 

K(2/) := C/„(y) - aia2?7„_i(y), n > 0. 

Hence, applying Theorem l3 . 31 and defining qo{x^ y) = 1, we see that the polynomials 

V„(y), A: = 0, 

qkix,y)Un~k(.y), l<k<n, 



Pn{x,y) 
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are orthogonal with respect to the weight function dfi{x, y) — dfiy{x)y^l — y^dy. 

Remark 5.1. These examples of Theorem [331 show that we can derive a complete 
basis of orthogonal polynomials in total order if N, the degree of h{z,y), is less 
than or equal to 4. There are other such weight functions, for example, those 
corresponding to 

h{z,y) = il-hz){l-b2z){l~2ayz + a^z^), |6i|,|62|<l, |a| < 1- 

In this case the integral d^y{x) yields a Bernstein-Szego weight with two linear 
factors in y (see Example 2). 
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